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Abstract 

We calculate the electric dipole moments of the neutron and the A 
within the framework of heavy baryon chiral perturbation theory. They 
are induced by strong CP-violating terms of the effective Lagrangian in the 
presence of the vacuum angle 6q. The construction of such a Lagrangian 
is outlined and we are able to give an estimate for 6q. 
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1 Introduction 



The axial U(l) anomaly in QCD implies an additional term in the Lagrangian, 
which violates P, T and CP. This new term is proportional to the so-called 
vacuum angle 8q, an unknown parameter, and its size may be determined from 
CP-violating effects, as e.g. 77 — > tttt or the electric dipole moments of the neutron 
and the A. The most recent measurements of the electric dipole moment of the 
neutron have constrained it to |l[ 

\dl\ < 6.3 x l(r 26 ecm. (1) 

On the other hand, theoretical estimates for <T n induced by the #o-term can De 
given leading to an upper bound for 9 Q 0, ||, [|] . Of particular interest here is the 
estimate of Pich and de Rafael who used an effective chiral Lagrangian ap- 
proach and came to the conclusion that it is possible to obtain an estimate for the 
size of the vacuum angle #0 with an experimental upper limit of \9q\ < 5 x 1CT 10 . 
However, the authors worked in a relativistic framework which does not have a 
systematic chiral counting scheme, so that higher loop diagrams contribute to 
lower chiral orders. This problem is avoided in heavy baryon chiral perturbation 
theory as proposed in 0], which allows for a consistent power counting. One 
should therefore study this effect within the heavy baryon formulation. Further- 
more, the baryon Lagrangian in |3| which describes the interactions of the neutron 
with the pseudoscalar nonet (ir, K, 77, 7/) does not contain explicitly CP-violating 
terms. They are rather induced by the vacuum alignment of the purely mesonic 
Lagrangian in the presence of the 8 -term. As shown in || |7j the most general 
baryonic Lagrangian taking the axial U(l) anomaly into account does have ex- 
plicitly CP- violating terms even at lower chiral orders. It has to be checked, if 
such terms lead to sizeable contributions for the present consideration. Finally, 
the authors of P| proposed to estimate the contribution from unknown counter- 
terms by varying the scale in the chiral logarithms. This procedure reveals the 
scale dependence of the involved coupling constants but not their absolute value, 
and it is desirable to have a somewhat more reliable estimate of the involved 
couplings. The aim of the present work is to reinvestigate the electric dipole 
moments of the neutron and the A by taking the above mentioned points into 
consideration. One has to check, if it is still possible to give a reliable estimate 
of the vacuum angle 6q, and if it is different from the one given in ||. 

The paper is organized as follows. In the next section we present the purely 
mesonic effective Lagrangian in the presence of the #o-term and the vacuum 
alignment is discussed. Baryon fields are included in the effective theory in Sec. 
3 by using the method outlined in ||. We proceed by calculating in Sec. 4 the 
electric dipole moment of the neutron and the A up to one-loop order within the 
heavy baryon framework. Numerical results and conclusions are given in Sec. 5. 
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2 The mesonic Lagrangian 



In this section, we will consider the purely mesonic Lagrangian in the presence 
of the #o-term. The derivation of this Lagrangian has been given elsewhere, see 



e.g. 0, |9], [HJ, so we will restrict ourselves to the repetition of some of the basic 
formulae which are needed in the present work. In |9|, the topological charge 
operator coupled to an external field is added to the QCD Lagrangian 

C = Cqcd ~ -£^9(x)tr c (G, u G^) (2) 

with G^ v = e^apG 01 ^ and tr c is the trace over the color indices. Under U(1)r x 
U(1) L the axial U{1) anomaly adds a term —(g 2 /16'K 2 )2NfatY c (G flu G flu ) to the 
QCD Lagrangian, with Nf being the number of different quark flavors and a the 
angle of the global axial U(l) rotation. The vacuum angle 9(x) is in this context 
treated as an external field that transforms under an axial U(l) rotation as 

9(x) -> 9'{x) = 9{x) - 2N f a. (3) 

Then the term generated by the anomaly in the fermion determinant is compen- 
sated by the shift in the 9 source and the Lagrangian from Eq. (|2]) remains invari- 
ant under axial U(l) transformations. The symmetry group SU(3)r x SU(3)l of 
the Lagrangian Cqcd is extendedQto U(3)r x U(3)l for C. The Green functions 
of QCD are obtained by expanding the generating functional around 9(x) = 9 
where the phase of the quark mass matrix emerging from the Yukawa couplings 
of the light quarks in the electroweak sector has been absorbed in 9q. The ex- 
tended symmetry remains at the level of an effective theory and the additional 
source 9 also shows up in the effective Lagrangian. Let us consider the purely 
mesonic effective theory first. The lowest lying pseudoscalar meson nonet is sum- 
marized in a matrix valued field U(x). The effective Lagrangian is formed with 
the fields U(x), derivatives thereof and also includes both the quark mass matrix 
M and the vacuum angle 9: C e Q(U,dU, . . . , M.,9). Under U(3)r x U(3)l the 
fields transform as follows: 

U' = RUL ] , M = RMI) , 9'{x) = 9{x) - 2N f a (4) 

with R G £7(3)^, L G U(3)l-, but the Lagrangian remains invariant. The 
phase of the determinant of U(x) transforms under axial U(l) as \ndetU'(x) = 
\ndetU(x) + 2iNf<y so that the combination 9 — ilndetU remains invariant. 
It is more convenient to replace the variable 9 by this invariant combination, 
£ e g(f/, dU, . . . , A4, 9—i In detU). One can now construct the effective Lagrangian 

3 To be more precise, the Lagrangian changes by a total derivative which gives rise to the 
Wess-Zumino term. We will neglect this contribution since the corresponding terms involve five 
or more meson fields which do not play any role for the discussions here. 
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in these fields that respects the symmetries of the underlying theory In partic- 
ular, the Lagrangian is invariant under U(3)r x £7(3) l rotations of U and M. at 
a fixed value of the last argument. The Lagrangian up to and including terms 
with two derivatives and one factor of M. reads 

+V i (UV fM U j )(U^U). (5) 

The expression (. . .) denotes the trace in flavor space and \ is proportional to 
the quark mass matrix \ = x — ZBqM. with M. = diag(m u , m^, m s ) and Bq = 
— (0\qq\0)/F% the order parameter of the spontaneous symmetry violation. The 
covariant derivative of U is defined by 

Vjj = dJj-i(v IM + a li )U + iU{v ll -ar). (6) 

The external fields v^x), a^{x) represent Hermitian 3x3 matrices in flavor space. 
Note that a term of the type iV 6 (WV ' Jj)V^6 can be transformed away [H and 
a term proportional to VJjV ll 9'W^9 does not enter the calculations performed in 
the present work and will be neglected. The coefficients Vi are functions of the 
variable 9 — zlndetC/, Vi(6 — ilndetU), and can be expanded in terms of this 
variable. The terms 14,.. .,4 are of second chiral order, whereas Vo is of zeroth chiral 
order. Parity conservation implies that the Vi are all even functions of 9—i In detC/ 
except V3, which is odd, and Vi(0) = V^O) = F%/4 gives the correct normalizaton 
for the quadratic terms of the Goldstone boson octet, where F n ~ 92.4 MeV is 
the pion decay constant. 

In order to use the effective Lagrangian, one must first determine the vacuum 
expectation value of U by minimizing the potential energy 

V(U) = V - V 2 (x ] U + W ] ) - iV 3 $U - xtf). (7) 

Since x is diagonal, one can assume the minimum Uq to be diagonal as well and 
of the form 

U = diag(e-^ , e~ iipd , e~ iips ) . (8) 
In terms of the angles <p q the potential becomes 

V(U ) = V (9o) - 4 B V 2 (9 ) £ m q cos cp g - 4 B V 3 (9 ) £ m o sin <Pv ( 9 ) 

9 9 

where we have introduced the notation 9 = 9 — J2 q fq- The Taylor expansions 
of the functions Vi read 

00 

r{2n)p2n 



Vtfo) = E Vi % n for i = 0,2 

n=0 

00 

v 3 (o ) = £K (2n+1) r +1 (10) 

n=0 
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with coefficients not fixed by chiral symmetry. Minimizing the potential with 
respect to the angles ip q leads to 



2B m q sin ip q = A + 2B B m q cos (p q 



with 



(JXJ 

A = 2fio(E^ (2n) ^o n 



n=0 



-1 



2B n 



-2nVi 2n) 6l n ~ l 



E 

n=l 

E m j c° s <Pj 

j=u,d,s 



-nv^ 2n) el n - 1 



-{2n-l)vt- 1] el n - 2 E rn jS m Vj ) 

j=u,d,s 



and 



OO -.OO 

B = ( E ^2 2n) #0") E ^3 2n+1) ^0 n+1 - 



n=0 



n=0 



To lowest order both in the quark masses m q and 1/N C Eq. QTTj) reads 
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-Bq Fl m q sin = Vq 6 
which is the equation for the tp q considered in j3], |TTJ. One then writes 

u = Ju Uy/u 

and U can be parametrized as 



(12) 



(13) 



(14) 



(15) 



U(<p,r)o) =exp{2i(p/F 7T + iJ-r] /F }, 



(16) 



where the singlet 7] couples to the singlet axial current with strength F Q . The 
unimodular part of the field U (x) contains the degrees of freedom of the Goldstone 
boson octet 6 



1 

72 



+ 



7T 



K 



K+ 
K° 



(17) 



75^ J 



while the phase detU(x) = e t ^ V0 ^ F ° describes the r] . The diagonal subgroup 
U(3)v of U(3)rxU(3)l does not have a dimension-nine irreducible representation 
and consequently does not exhibit a nonet symmetry. We have therefore used 
the different notation Fq for the decay constant of the singlet field. 
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One can now express the effective Lagrangian in terms of the Goldstone boson 
matrix U and the angles tp q 

= -Vo + V.iV^WU) + [V 2 + BV 3 ]( X (U + W)) -iAV 2 (U -rf) 
+ i[V 3 - BV 2 ]( X (U - rf)) + AV 3 (U + f/ f > + V,(UV^)(U^U),(18) 

where we have absorbed the angles (p q in hermitian matrices \ and H by defining 

Mx^o=X + iH, \Ju tfy/u~o = X-iH, (19) 

so that x — 2P diag(m g cos <p g ) and H = 2B diag(m q sin ip q ) = A + Bx- The 
Vi are functions of y/6r] /F + 9 Q , Vi(\/6r)o/Fo + 9 ) and we have assumed the 
external fields and a M to be diagonal which is the case if one considers elec- 
tromagnetic interactions. Note that the Goldstone boson masses at lowest chiral 
order are not only functions of the current quark masses m g , but also depend on 
the angles cos(p q . The kinetic energy of the 7/0 singlet field obtains contributions 
from VxCVpWVU) and V 4 (UV ^^(WV^U) which read 

(^ + ^^ 4 (0))^%. (20) 

We renormalize the rjQ field in such a way that the coefficient in brackets is 
1/2 in analogy to the kinetic term of the octet. By redefining F and keeping 
for simplicity the same notation both for r] and F one arrives at the same 
Lagrangian as in Eq. (0) but with V^O) = (F 2 — F%)/12 in order to ensure the 
usual normalization for the kinetic term of a pseudoscalar particle. 



3 CP- violating terms in the baryon Lagrangian 

As mentioned before, another source for CP-violation is the baryon Lagrangian. 
The CP- violating terms can be divided into two groups. Firstly, the vacuum 
alignment of the mesonic Lagrangian induces CP non-conserving meson-baryon 
interactions as considered in 0. But secondly, there are also explicitly CP- 
violating terms in the most general Lagrangian in the presence of the ^-vacuum 
angle which have been neglected in P|. In order to construct the Lagrangian in 
the baryon-sector, one has to adopt a non-linear representations for baryons. The 
main ingredient for a non-linear realization is the compensator field K(U, R, L) £ 
U(3)v, which appears in the chiral U(3)l X U(3)r transformation of the left and 
right coset representatives, £l(U) and £r(U): 

Il(U) - Li L (U)K\U,R,L) 

£ R (U) - Ri R (U)K\U,R,L). (21) 
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The field U from the last section is defined as 

u = inil 

The baryon octet B is given by the matrix 



(22) 



B 



V 



wo 



p 

n 



(23) 



which transforms as a matter field 

B -> B' = KBKl 
The covariant derivative of the baryon fields reads 

[D„B] = d,B+[T„B] 
with T M being the chiral connection 



(24) 



(25) 



(26) 



and r M = i^+a^, ^ 



-a M . For electromagnetic interactions the external fields 



are a M = and = —eQA^ with the quark charge matrix Q = |diag(2, —1, — 1). 
In order to incorporate the interactions with the mesons into the effective theory 
it is convenient to form an object of axial- vector type with one derivative 



iii = i [It - ir n) in - 3, - iQ II 
Further ingredients of the non-linear representation are 



t „~, 



(27) 



(28) 



and the quantity 

^ = &*2&± (29) 

where F^J L are the field strength tensors of r^/l^. The most general relativistic 
effective Lagrangian up to second order in the derivative expansion and contribut- 
ing to the electric dipole moments of the neutron and A reads 

C^b = iW 1 ([D^,B] ll ,B)-iW*(B ll ,[D^B}) + W 2 (BB) 

+W 3 {B^ 5 {e,B}) + W 4 {B^ l5 [e,B}) + W b (B llxlb B)(t) 
+tW 6 (B l5 B) + W 7 (B{x + , B}) + W 8 (B[x + , B}) + W 9 (BB) ( x+ ) 
+iW w (B{x-, B}) + iW 11 (B[x-,B}) + iW 12 (BB)(x_) 
+iW 13 (Ba^ 5 {F£ u , B}) + iW 14 {Ba^ 5 [F^, B]) 



+iW 15 {Ba^ 5 B){F: 



n 



(30) 
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The Wi are functions of the combination y/Erjo/Fo + 0$. From parity it follows 
that W\ t ...$ and W 7i8) 9 are even in this variable, whereas Wq and Wio,...,i5 are 
odd. The latter have not been taken into account in j3]. One can further reduce 
the number of independent terms by making the following transformation. By 
decomposing the baryon fields into their left- and right handed components 

B R/L = ^(l± l5 )B (31) 

and transforming the left- and right-handed states separately via 

Br/l ~* VW^W 6 Br/l 

Br/l Vw 2T iW 6 BR/L (32) 

one can eliminate the (B^B) term and simplify the coefficient of (BB). The 
details of this calculation are given in || . Note that this transformation leads to 
mixing of the terms of the type B\±B with B^\±B which are of third chiral 
order and have been neglected here. Furthermore, the terms Wn,...,i& mrx with 
the terms from the Lagrangian 

C = W w (Ba^{F^, B}) + W 17 (Ba^[F^, B]) + W 18 (Ba^B)(F^), (33) 

but in both cases the form of the Lagrangian does not change and we can proceed 

o o 

by neglecting the W% term and setting W% = W* and W2 = — M with M being 
the baryon mass in the chiral limit . The expansion of the coefficients in terms 
of the Wi read 

W 1 = -- + ... , W 2 = - M, 

W 3 = ~D + ..., w A = ~F + ..., W 6 = ±\ + ..., 
W 7 = b D + ..., W 8 = b F + ..., W 9 = b + ... } 

W l = w l (^ rVo + e ) + ..., for z = 10,...,15 , (34) 

-TO 

where the ellipses denote higher orders in \ZQt]q /Fq + % and we have only shown 
terms that can contribute to one-loop order. The axial- vector couplings D and F 
can be determined from semileptonic hyperon decays. A fit to the experimental 
data delivers D = 0.80 ± 0.01 and F = 0.46 ± 0.01 0. The third coupling, A, is 
specific to the axial flavor-singlet baryonic current. The coefficients bo,F,o have 
been determined from the calculation of the baryon masses and the ttN a-term 



up to fourth chiral order [13| and their mean values are, in units of GeV 



>D 



0.079, b F = -0.316, b = -0.606. (35) 
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The numerical values of the parameters u>io,...,i5 are not known. We leave them 
undetermined for the time being and will give later an upper bound. 
Again the correct vacuum has to be chosen. This is done by setting 



R 

and we will choose the coset representatives such that 

t R = & = u = y/U. 
For diagonal external fields and a M one can write 



and 



r = r 

1 H t 1 



u 



t (dp - ir^) u + u(d tl - il^) u ] 



u ] (<9 M - ir^) u - u(d^- il/j) u ] 



Furthermore, one obtains 

X- = X--iA(U + tf)-iB X+ , 
where the quark mass matrix enters in the combinations 

X± = ux ] u±u ] xu ] - 

Finally, simplifies to 

K = F tv = + «F>t. 

The relativistic baryon Lagrangian reads to the order we are working 



(36) 
(37) 

(38) 
(39) 

(40) 
(41) 
(42) 



C<PB = i(B^[D»,B])-M(BB) 

~D(B W {u»,B}) - ±F{B W [u»,B}) + ±\(B W B)(u») 

-ib D A{B{U - U\ B}} - ib F A{B[U - U\ B\) - ib A{BB)(U - U*) 

/fi /fi 
+AAw w ^rVo(BB) + 6Aw 12 ^rVo(BB) 

-TO 

+i(w' 13 9 + w^^iBa^iF^, B}) 

-TO 



+i(w[ 4 9 + w u — rj ) (Ba^ [F^, B]) 

-TO 

+i(w[ 5 9 + w 15 ^-r] )(Ba^j 5 B) (F^) , 



(43) 
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where we have neglected meson-baryon interactions with more than one meson 
field since they do not contribute at one- loop order and terms of O(0q) or higher 
orders are omitted throughout this work. Note that there exist terms of fourth 
chiral order of the type Ba llv ^X-F+ B. Using Eq. (fHf) for \- they induce CP- 
violating terms of the form OqBo ^75-^+^ ' B which are already accounted for by 
the terms 1013,14,15. This amounts to a renormalization of the couplings ^0^13,14,15- 
We have therefore introduced the notation w' 13 14 15 for these interaction terms, 
in order to distinguish them from the unrenormalized 1013,14,15 of the interactions 
proportional to r] . 

The drawback of the relativistic framework including baryons is that due to 

o 

the existence of a new mass scale, namely the baryon mass in the chiral limit M, 
there exists no strict chiral counting scheme, i.e. a one-to-one correspondence 
between the meson loops and the chiral expansion. In order to overcome this 
problem one integrates out the heavy degrees of freedom of the baryons, similar 
to a Foldy-Wouthuysen transformation, so that a chiral counting scheme emerges. 
To this end, one constructs eigenstates of the velocity projection operator P v = 
(l + V)/2 

B v (x) = e^ IV - x P v B(x). (44) 

The Dirac algebra simplifies considerably. It allows to express any Dirac bilinear 
B V T^B V (r M = 1,7^,75,...) in terms of the velocity and the spin operator 
2S^ = i^a ' lw v v . One can rewrite the Dirac bilinears which appear in the present 
calculation as 

B v j^ 5 B v = 2B V S^B V , B v a^ u j 5 B v = 2i (v fl B v S u B v - v u B v S^B v ). (45) 

In the following, we will drop the index v. The Lagrangian of the heavy baryon 
formulation reads 

C+b = i{B[v ■ D,B]) - D{BS^{u^B}) - F(BS^,B]) + X(BS fl B)(u' i ) 

-ib D A{B{U - XJ\ B}} - ib F A{B[U - U\ B]) - tb A{BB}(U - U ] ) 

/fi /fi 
+AAw 10 ^ rVo (BB) + 6Aw 12 ^ r r ]o (BB) 

-TO -TO 

\/6 

-4K 3 O + wu^-TtiiBvpSviFP, B}) 
-ro 

\/6 

-4«0 O + w^^iBv^iF^, B]) 
^0 

-A(w' 15 6 + wV^Bv^B)^). (46) 
^0 

No relativistic corrections are needed to the order we are working. 
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4 The electric dipole moments and 

In this section, we will calculate the electric dipole moments of the neutron and 
the A at lowest order in chiral perturbation theory, i.e. 0(p 2 ). Due to the vacuum 
alignment the baryon Lagrangian contains interaction terms &d,f,o of zeroth chiral 
order and one-loop diagrams with these vertices will also contribute at 0(p 2 ). In 
the relativistic framework of the electric dipole moment of the neutron d? n has 
been defined via ^ 

£nEDM = -dlnia^nF^, (47) 

where F^ u is the field strength tensor of the photon field A^. We prefer to rewrite 
this form factor 

DZ(q 2 )u(p'W»lMpW (48) 

with q = p' — p being the momentum transfer. The electric dipole moment is 
given by 

dr = Dytf = 0). (49) 

For the calculation of the form factor in the heavy baryon approach we set v M = 
(1,0) and use the Breit frame v ■ p = v ■ p' since it allows a unique translation 
of Lorentz-covariant matrix elements into non-relativistic ones. In this frame the 
form factor reads 

- 2iDl(q 2 )Hv u S ■ qH + . . . , (50) 

where H is the large component of u and the ellipsis stands for a similar expression 
in the small components of u, which is of higher chiral order and will be omitted. 
The electric dipole moment d% receives contributions from the u>i3-term and the 
loops (we work in the isospin limit m u = = m) 

dl = dl {tree) + dl {loov) (51) 

with 

d ^ree) = _g g ^ [ 1 ^ + —^—V^V^ } , (52) 

7T 7]0 

where the pertinent diagrams are shown in Figure 1. Diagram lb) is missing 
in H since the mesonic Lagrangian used within this work does not have a term 
linear in the singlet field 7] . Such an interaction originates from the term V^(x-) 
which is not considered in 0. The chiral logarithms of the diagrams shown in 
Fig. 2 read 

dl^ = ^eV i2) ^-(b D +b F )(D+F)ln^ + (b D -b F )(D~F)ln^f] (53) 



with fi being the scale introduced in dimensional regularization. In the Breit 
frame only diagrams 2a) and b) contribute to the electric dipole moment. Di- 
agrams 2c), d) vanish and 2e), f) are proportional to S v and therefore do not 
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contribute to d%- The loop integral for diagrams 2a) and b) contains also ana- 
lytic and divergent pieces which can be absorbed by redefining w' 13 . The divergent 
pieces of w' 13 cancel the divergencies from the loops and render the final expres- 
sion finite. We summarize the remaining analytical contributions in w{ 3 , so that 



w[ 3 in Eq. 



is understood to be finite. The results in Eqs. 



and (p31) are 



not in contradiction with the fact that the electric dipole moment induced by the 
#o-term tends to zero if any of the quark masses vanish. If, e.g., m u = then a 
solution for Eq. (|TTD is given by ip u = 9 and <fd,s — leading to 9 = 9 — ip u = 0. 
Therefore, vanishes in this case. 
The results for the A are 



d 



7 (tree) 
A 



-4e9 r 



1 r 



+ 



16 



F 2 Fnm 2 

1 TT 1 " l r)0 



1/ (1 V (2) 7„ 



and 



d 



7 (loop) 
A 



(2) 



9 [b D F + b F D] In 



(54) 



(55) 



Note that the relation d^ n = 2d\ is only valid at tree level and not for the chiral 
logarithms as claimed in ||. The discrepancy is due to the lack of pion loops in 
the present work. Once one accounts for the mistake made in [[| by replacing 
In by In mx in the chiral loop contribution for d\ one obtains our result 



5 Numerical results and conclusions 

In order to compute the numerical results for the electric dipole moments shown 
in the last section, we use the central values for the parameters bo = 0.079 GeV -1 , 
bp = —0.316 GeV -1 , D = 0.80 and F = 0.46. To lowest order in the angles ip q 
and using rh -C m s one can express 9q in terms of 9q via 

OT/(2) 

9 ~ [1 + (56) 

F 2 m 2 

From the calculation of the r] and rj masses and decay constants one can extract 
the value for [[TJ 

V {2) ~ -^F 4 ~ -5.0 x 10- 4 GeV 4 , (57) 



so that 

j?V 

9 ^-^f ^o^-0.04^o. (58) 

Inserting this into the loop contribution and using fi = 1 GeV, m m ~ m v > = 958 
MeV we obtain 



d i(ioo P ) = _7. 5 x io- 16 ^ ecm 

d\ {loop) = -1.7 x Kn 16 o ecm. (59) 
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The numerical result for d^ loop ^ is in agreement with the one given in || once 
one accounts for the different values of the parameters bp, bp, D, F and the scale 
fj, used within that work. The result for the chiral logarithm of d\ ^ is consid- 
erably smaller than for d^ loop ^ since there is no contribution from the pion loops 
which dominate in the case of d%. 

A precise numerical value for the tree contribution to the electric dipole mo- 
ments cannot be given since the parameters 1013 and w[ 3 are not known. However, 
we will give an upper bound for their contribution based on large N c arguments 



which seem to work well in the purely mesonic sector |14[ [Lql . In the present 
investigation we are only interested in an order of magnitude estimate for the 
vacuum angle 9 and for this purpose it is sufficient to give a numerical range for 
the tree level contribution. We will first estimate the ratio of diagrams la) and 
lb) using l/N c arguments. Applying large N c counting rules, see e.g. [|, [10[| , 
both u> 13 an d w[ 3 are of order 0(N®) so that we can assume |xf 13 / 1^13 1 = Oil). 
One obtains the ratio 

« %(V a, N(U2i (60) 



\ r n( la )\ u^ la )| F 4 m 2 , 

\Un I I it rf 

where we have used F n /F = 1 + 0{N^ 1 ) and taken the value for V3 from |14 



V 3 {1) ~ 0.04F* ~ 3.5 x 10" 4 GeV 2 . (61) 

Diagram lb) turns out to be insignificant in our estimate. 

Based on large N c arguments one can also give an upper bound for w[ 3 . The 
contact terms of the Lagrangian in Eq. (|33|), which contribute to the magnetic 
moments of the baryons, describe - similar to the w[ 3 -term - the coupling of the 
field strength tensor F^ v to the baryons. But the leading coefficient of the Taylor 
expansion of Wiq, Wiq, is of order O(N^), whereas \w[ 3 \ ~ \wi 3 \ = 0(N°), so that 
we can assume \wi 3 \ < \wiq\. In a calculation of the baryon magnetic moments 
to fourth chiral order wiq has been determined to be |16| 



w 16 ~ 0.4 GeV . (62) 
Inserting this upper limit for w[ 3 one obtains 

\d? n {tree) \ < 9.6 x 10- 16 # ecm 

\d\ {tree) \ < 4.8 x 10- 16 ^ ecm. (63) 

Since we have taken quite conservative limits, the tree level contributions could 
be dominating if the extreme values are chosen. A more realistic estimate might 
be obtained by setting \wi 3 /w\q\ ~ | for N c = 3, so that 

(*r«0 1 ~ 3 . 2 X HT 16 O e Cm 

\d\ {tree) \ ~ 1.6 x 10- 16 ^ ecm. (64) 
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While the chiral logarithms dominate for d^, the tree level contribution can be 
substantial for the A. This leads to 

dl = (-7.5 ± 3.2) x KT 16 9 e cm 

(f A = (-1.7 ± 1.6) x 10~ 16 # ecm, (65) 

where the theoretical uncertainty is given by the estimate of the tree level con- 
tribution in Eq. (0). From a comparison of the central value for (P n with the 
experimental upper limit in Eq. (|l|) we obtain 

|0 O | < 8.4 x KT 11 . (66) 

In the case of the A the experimental constraint is given by 

(f A < 1.5 x KT 16 ecm (67) 

which leads to 

|0 O | < 0.9. (68) 

We have shown that it is possible to obtain a reliable limit for the vacuum angle 
8q by calculating the electric dipole moment of the neutron within the framework 
of heavy baryon chiral perturbation theory. To this end, we have constructed 
the most general effective Lagrangian up to one-loop order in the presence of the 
vacuum angle 8q with the method proposed in ||. The theoretical uncertainty 
from unknown parameters at tree level has been estimated by using large N c 
arguments. While the chiral loops dominate for the electric dipole moment of the 
neutron, the counterterms can be substantial in the case of the A. 
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Figure captions 



Fig.l Shown are the tree diagrams for the electric dipole moment. Solid and 
dashed lines denote baryons and pseudoscalar mesons, respectively. The 
wavy line represents a photon and the dot is a CP- violating vertex. 

Fig. 2 Loop diagrams contributing to the electric dipole moment. Solid and dashed 
lines denote baryons and pseudoscalar mesons, respectively. The wavy line 
represents a photon and the dot is a CP- violating vertex. 
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Figure 1 






Figure 2 
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